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We study the extent of the conformal window for an SU(3) gauge theory with Nt Dirac fermions in the 
fundamental representation. We present lattice evidence for 12 < Nf < 16 that the infrared behavior is 
governed by a fixed point, while confinement and chiral symmetry breaking are present for Nf < 8. 
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With a small number of massless fermions, a vector-like 
gauge field theory such as QCD exhibits confinement and 
dynamical chiral symmetry breaking. But if the number of 
massless fermions, Nf, is larger, near but just below the 
value, Nf f , at which asymptotic freedom sets in, the theory 
is conformal in the infrared, governed by a weak infrared 
fixed point (IRFP) which appears already in the two-loop 
beta function lily]. There is no confinement, and chiral 
symmetry is unbroken. It is thought that this IRFP persists 
down to some critical value Nf, where the coupling is suf- 
ficiently strong that the transition to the confined, chirally 
broken phase takes place. The range A^ f > Nf > Nj is 
the "conformal window", where the theory is in the "non- 
Abelian Coulomb phase". 

Theories in or near the conformal window could play a 
key role in physics beyond the standard model. For ex- 
ample, a theory near the conformal window could describe 
electroweak symmetry breaking. It is therefore important 
to study the extent of this window, as well as the order of 
the transition at Nj and the properties of the theory inside 
the window and near it. Despite interest in these questions 
for many years, little is known with confidence. This can 
be contrasted with supersymmetric QCD, where duality ar- 
guments determine the extent of the conformal window 
and lead to weakly coupled effective low-energy theories 
at both ends H. 

An upper limit on Nf for both supersymmetric and non- 
supersymmetric theories has been proposed based on the 
counting of massless degrees of freedom, employing the 
thermodynamic free energy [4]. For a supersymmetric 
SU(A) gauge theory with Nt massless Dirac fermions in 
the fundamental representation (where Nf = 3 AO, one 
finds Nf < (3/2) A, a limit precisely saturated by the 
result from duality arguments. For a non-supersymmetric 
SU(A) theory with Nf massless Dirac fermions in the 
fundamental representation (where Aj! f = (11/2) A), one 
finds Nf < 4A[1 - (1/18A 2 ) + ...]. It is not known to 
what extent this limit is saturated. The most recent lattice 
studies 15|], for A = 3, lead the authors to the conclusion 
that Nf is much lower. They find 6 < A^ < 7. 

In this letter, we describe a new lattice study of the con- 
formal window for an SU(3) gauge theory with Nf Dirac 
fermions in the fundamental representation. We adopt a 
gauge-invariant definition of the running coupling, valid 



for any strength, derived from the Schrodinger functional 
(SF) of the gauge theory 0, 0]. For an asymptotically 
free theory, this coupling agrees with the perturbative run- 
ning coupling at short enough distances Hi. 

Making use of staggered fermions as in Ref. OIL 

it is 

most straightforward to restrict attention to values of A^ 
that are multiples of 4. The value Nf = 16 leads to an 
IRFP that is sufficiently weak that it is best studied in per- 
turbation theory. The value A^ = 4 is expected to be well 
outside the conformal window, leading to confinement and 
chiral symmetry breaking as with Nf = 2. Thus we focus 
on N f = 8 and N f = 12. 



The Schrodinger functional is the transition amplitude 
from a prescribed state at time t = to another state at 
time t = T. It can be written as a Euclidean path integral 
in a spatial box of size L with Dirichlet boundary condi- 
tions at t = and t = T where T is 0(L). Periodic 
boundary conditions are imposed in the spatial direction. 
The Schrodinger functional can be written as 
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where U are the gauge fields and X are the staggered 
fermion fields. W and W are the (fixed) boundary values 
of the gauge fields, and £, £ are the boundary values 
of the fermion fields at t = and t = T, taken here to be 
zero. The quantity Sq is the Wilson gauge action and Sp 
is the massless staggered fermion action. 

The gauge boundary values W(rj), W (77) are chosen 
such that the minimum action configuration is a con- 
stant chromo-electric field fill [l2ll whose magnitude is of 
Q(l I V) and is controlled by a dimensionless parameter r/ 
[fl3ll . The SF running coupling g 2 (L, T) is defined by tak- 
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where k = 12 (£) [sin (2vra 2 /3LT) + sin (vra 2 /3LT)] 
is chosen so that (7 2 (L, T) equals the bare coupling at tree 
level. In general, g 2 (L,T) measures the response of the 
system to small changes in the background chromo-electric 
field. 

For staggered fermions, L /a must be even but T/a must 
be odd, where a is the lattice spacing. To cancel the resul- 
tant 0(a) bulk lattice artifact, the coupling is defined as the 



2 



average over T = L ± a: 
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g 2 (L,L-a) f(L,L + a 
0(a) terms on the Dirichlet boundaries remain 

01. We 

include a perturbative one-loop counterterm of O(g^a) in 
our calculations to remove partially the O(o) boundary ar- 
tifact. Since g 2 (L) depends on only one IR scale, L, it 
provides a technical advantage for studying an IRFP over 
other possible non-perturbative definitions of the running 
coupling, such as from the static potential V(r), which 
must be computed from Wilson loops at scales r <C L 
to avoid finite-size effects [fl4fl . 

To set the stage, we review briefly the behavior of g 2 (L) 
in continuum perturbation theory through three loops. By 
computing ~g 2 (L) in lattice perturbation theory and setting 
to zero terms that vanish as a/L — > 0, a continuum beta 
function, j3, may be defined such that L(d / dL)~g 2 (L) = 
P [f(L)) = b x t(L) + b 2 t(L) + 6 3 7J 8 (L) + • • • , where 
the first two (scheme-independent) coefficients are 



(4tt) 



[11 



iNf] 



(4*+ 



[102 - f N f ) . (4) 



The third coefficient is scheme dependent, given in the SF 
scheme by |@] 
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+ 1.197(10) + 



6 J 2vr 8vr 2 
where c 2 = 1.256 + 0.040^ andc 3 = 6 
0.140(6)^ - 0.0330(2)^), and where feipis the three- 
loop coefficient defined in the MS scheme (only in the loop 
expansion), given by 
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For N-f = 16, a weak two-loop IRFP exists at g t ~ 
0.52. The higher order corrections are very small in both 
the SF and MS schemes. 

For Nf = 12, the two- and three-loop beta functions 
also exhibit an IRFR although the reliability of the loop 
expansion is less clear. The two-loop IRFP is at g 2 ~ 
9.48. At three loops the fixed point strength is reduced by 
roughly 50%, to g\ ~ 5.18 (5.47) in the SF (MS) scheme. 
(In the MS scheme, where a four-loop result exists, the 
fixed point value increases slightly to 5.91.) Since the cor- 
responding loop expansion parameter gl/Air 2 is smaller 
than unity, perturbation theory could provide a reasonable 
basis of comparison for our lattice simulations. 

For N f = 8, there is no two-loop IRFP. While an IRFP 
can appear at three loops and beyond, its scheme depen- 
dence and typically large value means that there is no ev- 
idence for an IRFP accessible in perturbation theory. A 
non-perturbative study is essential. 

With this as background, we next describe our lattice 
simulations for Nf = 8 and 12. For Nf = 8 we compute 
g 2 (L) for bare lattice couplings = 6/g 2 £ [4.5,7.1] 
with typically 0.1 spacing. For Nf = 12 we choose 



[3 G [4.15,6.5] also with typically 0.1 spacing 12811 . In 
both cases, we use lattice extents L/a = 4, 6, 8, 10, 12, 16 
and 20, with the larger L/a computations done at fewer 
values due to the much higher computational cost. We use 
the standard hybrid molecular dynamics (HMD) R algo- 
rithm [15] with unit length trajectories. We generate three 
independent ensembles varying the number of steps per tra- 
jectory, typically in the range 64-128 steps but up to 512 
steps at stronger couplings, and perform a quadratic extrap- 
olation to remove finite step size errors. Within our statis- 
tics, the observed systematic shift in ~g 2 (L) due to finite 
step size is negligible over the chosen range of step sizes. 

For the majority of the simulations, those at relatively 
weak lattice coupling, we employ of order 40,000 HMD 
trajectories, sufficient to estimate reliably autocorrelations. 
At stronger couplings, the autocorrelations become much 
longer, with the time histories showing the previously ob- 
served phenomenon 1 13 , of large excursions lasting a 
few thousand trajectories. In such cases, simulations are 
run longer, up to 80,000 HMD trajectories. Error estima- 
tion is performed using the jac knife method with the block 
size adjusted to eliminate the effects of autocorrelations. 

To observe the running of the coupling over a large range 
of scales requires the generation of several hundred inde- 
pendent ensembles at various values of the box size, L/a, 
and bare gauge coupling, [3 = 6/ 1 g 2 . With so many in- 
dependent statistical estimates of ~g 2 (L), occasional large 
statistical fluctuations of these estimates are expected. So, 
we model our estimates with a smooth interpolating func- 
tion based on a truncated Laurent series 



g 2 ([3,a/L) = Y, 
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UP-Po(a/L)] K 
with polynomial dependence of the coefficients on a/L 

Best-fit values for the coefficients at 8 and 12 flavors will 
be included in a future paper. This interpolating function 
is used only to describe the lattice data in the limited range 
where the data exists, well away from the poles. 

The extrapolation to the continuum is implemented 
using the step-scaling procedure |17, OH], a systematic 
method that captures the renormalization group evolution 
of the coupling in the continuum limit. The basic idea is 
to match lattice calculations at different values of a/L, by 
tuning the lattice coupling = 6/ ' g 2 so that the coupling 
strength g 2 (L) is equal on each lattice. Keeping ~g 2 (L) 
fixed while changing a/L allows one to extract the lat- 
tice artifacts. Previous work by the ALPHA collabora- 
tion has shown that perturbative counterterms greatly re- 
duce 0(a/L) artifacts from the running coupling. So, our 
preferred method of extrapolating to the continuum limit at 
each step assumes that 0(a 2 /L 2 ) errors dominate. 

In practice, one calculates a discretized version of the 
running coupling, known as the step-scaling function. In 
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the continuum, it is designated 

a(s,g 2 (L)) = g 2 (sL), (9) 

where s is the step size. On the lattice, a/L terms are also 
present; one defines S(s, ~g 2 (L), a/L) similarly, so that it 
reduces to a(s,g 2 (L)) in the continuum limit: 

a(s,g 2 {L)) = lim£( S ,g 2 (L),a/L) (10) 

a— >0 

First, a value u = g 2 (L) is chosen. Several ensembles with 
different values of a/L are then generated, with j3 = 6/g 2 
tuned using our interpolating function Eq. (0) so that the 
measured value of g 2 (L) = u on each. Then for each (3, 
a second ensemble is generated with larger spatial extent, 
L — > sL. The measured value of g 2 (sL) on the larger 
lattice is exactly X(s,u, a/L). Extrapolation in a/L to 
the continuum then gives us the value of a(s, u). Using 
the value of a(s, u) as the new starting value, this process 
may then be repeated indefinitely, until we have a series of 
continuum running couplings with their scale ranging from 
L to s n L. In this paper, we take s = 2. 

Our results for Nf = 12 continuum running are pre- 
sented in Fig. [TJ We take L to be the scale at which 
~g 2 (L) ~ 1.6, a relatively weak coupling. The points 
shown are for values of L/L increasing by factors of 
two. The step-scaling procedure leading to these points 
involves stepping L/a from 4 — > 8, 6 — > 12, 8 — > 16, 
and 10 —* 20, and then extrapolating S(2, u, a/L) to the 
continuum limit assuming that 0(a 2 /L 2 ) terms dominate. 

Various sources of systematic error must be accounted 
for. The interpolating function Eq. (0) may not contain 
enough terms to capture the true form of ~g 2 {L) at large 
L where there is sparse data, and although the 0(a/L) 
terms are expected to be small, ignoring them completely 
in the continuum extrapolation may introduce a small sys- 
tematic effect. In addition, a few simulations that had run at 
least 20% of their target length, but were not yet completed, 
were included in the fit. The statistical error in g 2 (L) for 
these cases was likely underestimated. Here we provide 
an estimate of our systematic error by varying our con- 
tinuum extrapolation method between extremes. Inspec- 
tion of T*(2,u,a/L) as a function of a/L indicates that 
dropping the step from 4 — > 8 and performing a constant 
extrapolation underestimates the true continuum running, 
while performing a linear fit to all four steps gives an over- 
estimate. These define the upper and lower bounds of the 
shaded region in Fig.Q] which we take to be a conservative 
estimate of the overall systematic error. 

The observed IRFP for Nf = 12 agrees within the 
estimated systematic-error band with three -loop perturba- 
tion theory in the SF scheme. An important feature for 
Nf = 12 is that the interpolating curves are anchored by 
values of g 2 (L) that are also above the IRFP. For < 4.4, 
~g 2 (L) is large, decreasing as L/a increases with fixed j3. 
In the step-scaling function, values of u in this range lead 
to S(2, u, a/L) < was a/L — > 0. This behavior is similar 
to that found in Ref. [19] for Nf = 16, and consistent with 
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FIG. 1: Continuum running coupling from step scaling for Nf = 
12. The statistical error on each point is smaller than the size of 
the symbol. Systematic error is shown in the shaded band. 



approaching the IRFP from above in the continuum limit. 
In a future paper, we will exhibit both the step-scaling re- 
sults and the continuum evolution in this region. 

Our results for Nf = 8 continuum running are presented 
in Fig. |2j starting at a scale L where ~g 2 (L) ~ 1.6, and 
exhibiting points with statistical error bars for values of 
L/Lq increasing by factors of two. The three step-scaling 
procedures are the same as in the Nf = 12 case. Step- 
ping L/a from 4 -» 8, 6 -> 12, 8 -> 16 and 10 -> 20 
with quadratic extrapolation again provides the points with 
statistical error bars, and the other two procedures define 
the upper and lower bounds of the systematic-error band. 
For comparison, we have also shown the results of two- 
and three-loop perturbation theory up to ~g 2 ~ 10, beyond 
which there is no reason to trust the perturbative expansion. 

The Nf = 8 running coupling shows no evidence for an 
IRFP, or even an inflection point, up through values exceed- 
ing 14. The points with statistical errors begin to increase 
above three-loop perturbation theory well before this value. 
This behavior is similar to that found for the quenched the- 
ory Q and for N f = 2 iH, although, as expected, the 
rate of increase is slower than in either of these cases. The 
coupling strength reached for Nf = 8 exceeds rough es- 
timates of the stren gth requ ired to trigger dynamical chiral 
symmetry breaking l2lll22L 12311 . and therefore also confine- 
ment. This conclusion must be confirmed by simulations 
of physical quantities such as the quark-antiquark potential 
and the chiral condensate at zero temperature. 

To conclude, we have provided evidence from lattice 
simulations that for an SU(3) gauge theory with Nf 
Dirac fermions in the fundamental representation, the value 
Nf=8 lies outside the conformal window, and therefore 
leads to confinement and chiral symmetry breaking; while 
Nf =12 lies within the conformal window, governed by an 
IRFP. We stress that these conclusions do not depend cru- 
cially on the L/a = 4 data, which are of limited use in 
the SF scheme |gj. Thus the lower end of the conformal 
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FIG. 2: Continuum ranning coupling from step scaling for Nf — 
8. Errors are shown as in Fig.Q] Perturbation theory is shown up 
to only g 2 (L) ~ 10. 



window, Nf, lies in the range 8 < N c f < 12. 



This conclusion, in disagreement with Ref. |5|], is 
reached employing the Schrodinger functional (SF) run- 
ning coupling, ~g 2 (L). This coupling is defined at the box 
boundary L with a set of special boundary conditions. It 
runs in accordance with perturbation theory at short enough 
distances, and is a gauge-invariant quantity that can be 
used to search for conformal behavior, either perturbative 
or non-perturbative, in the large L limit. 

For Nf=8, we have simulated g 2 (L) up to values that 
exceed rough estimates of the coupling strength required 
to trigger dynamical chiral symmetry breaking |2U l22l 



23l 12411 . with no evidence for an IRFP. For Nf=l2, our 
observed IRFP is rather weak, agreeing within the esti- 
mated errors with three-loop perturbation theory in the SF 
scheme. 

The simulations of g 2 (L) at Nf=8 and 12 should be 
continued to achieve more precision. It is also important 
to supplement the study of ~g 2 (L) by examining physical 
quantities such as the static potential, and demonstrating 
directly that chiral symmetry is spontaneously broken for 
Nf=8 through a zero-temperature lattice simulation. Sim- 
ulations of ~g 2 (L) for other values of Nf, in particular 
Nf=10, will be crucial to determine more accurately the 
lower end of the conformal window and to study the phase 
transition as a function of Nf. All of these analyses should 
be extended to other gauge groups and other representation 
assignments for the fermions [25]. 
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